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The step edge barrier of a prototypical organic semiconductor molecule, 3,4,9, 10-perylene- 
tetracaboxylic-dianhydride (PTCDA) has been analysed by means of calculations based on em- 
perical potentials. The minimum energy path (MEP) has been calculated for a single molecule 
on a substrate of three molecular layers between equivalent minimum energy positions within two 
neighboring unit cells. To determine the step edge barrier, we have calculated the MEP over a step 
to a fourth layer of molecules. We found energy barriers of Ed = 80 meV for in-layer diffusion and 
Es = 750 meV for step crossing, indicating a strong Ehrlich-Schwoebel effect for PTCDA. 

PACS numbers: 81.15.Aa, 68.43.Jk 61.66. Hq 
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In the last decades, there has been a growing inter- 
est in electronic devices built with organic semiconduc- 
tor molecules. Organic molecular devices, in contrast to 
present-day electronics based on silicon, offer a higher 
mechanical flexibility and lower production costs. How- 
ever, this field of research still faces many open questions 
and future challenges. One of them is the understanding 
of the growth of organic crystals^ as the quality of an or- 
ganic crystal has a strong influence on device properties 
such as, e.g., the charge carrier mobility. 

An important factor in inorganic crystal growth ki- 
netics is the Ehrlich-Schwoebel (ES) effect, 2 '? which 
refers to the additional energy barrier experienced by an 
atom or molecule crossing a step edge. The strength 
of the step edge barrier controls the amount of mass 
transport between different crystalline layers, and hence 
determines the temperature dependence of the kinetic 
growth mode (two-dimensional layer-by-layer or three- 
dimensional mound growth) in situations where energetic 
driving forces shaping the surface morphology can be 
neglected A 

There is considerable indirect evidence that points 
to the presence of a step edge barrier for prototypical 
organic semiconductor molecules like 3,4,9, 10-perylene- 
tetracaboxylic-dianhydride (PTCDA)£&£££ and 
pentacene ) 10 ' 11 ' 12 and the ES effect has been invoked in 
several simulation studies of organic film growth ^ 10 ' 13 . 
Very recently a quantitative estimate for the step edge 
barrier on pentacene has been obtained from an analysis 
of experimental layer coverages^. So far, however, the 
existence of step edge barriers for organic molecules has 
not been established on the molecular level. The problem 
in theoretical studies is the high computational effort: 
Due to the size of the relevant systems («50 atoms 




FIG. 1: PTCDA molecule (C240 6 H 8 ). The size of the 
molecule (with van der Waals radii) is 14.2 Ax 9.2 A. The 
two oxygen atoms denoted by OS have slightly different pa- 
rameters in the AMBER force field than those denoted by 
O. 



per molecule for relevant molecules, «100 molecules in 
a model molecular crystal), ab initio methods are not 
suitable. Empirical potentials offer an opportunity to 
avoid these problems, if one accepts the corresponding 
loss in accuracy. The capability of empirical poten- 
tials for growth studies of organic thin films has been 
demonstrated in recent publications ! 15 ' 16 ' 

In this paper, we present calculations of the diffusion 
and step edge barriers for a PTCDA molecule (see Fig. 
[lj on a PTCDA crystal. Energies have been calculated 
using empirical potentials; the general method presented 
here should also apply to other molecular systems. 

To model a PTCDA single crystal surface, a rigid cell of 
three layers of PTCDA molecules in a-configuration was 
set up, using the crystallographic data by Ogawa et alJ£. 
The optimum postion for a single rigid molecule on the 



topmost layer was determined using a gradient search al- 
gorithm. The diffusion barrier Ejy was then calculated by 
finding the minimum energy path (MEP) between mini- 
mum positions within two neighboring PTCDA unit cells, 
using the nudged elastic band algorithm^ Potential en- 
ergies have been calculated by means of empirical pair 
potentials of Lennard-Jones (LJ) type. Parameters from 
the AMBER molecular force field^ were used. AMBER 
uses a 6-12 LJ potential with the energy given by 
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where i and j denote all pairs of atoms with distance 
Rij in the PTCDA substrate (i) and the single PTCDA 
molecule on top (j). Aij and Bij are specified by the 
parameters r* and e which are given for each atom type 
(see table U). The combination rules are 
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PTCDA is a polar molecule, therefore we had to ac- 
count for electrostatic interactions. The electrostatic en- 
ergy has been calculated by assigning partial charges to 
all atoms by a Mulliken method, using the commercial 
software Hyperchem. A dielectric constant of e r = 3.2 
was used, which was found as an average value for bulk 
PTCDA. 21 The total electrostatic energy is then calcu- 
lated by summing up all pairs of charged atoms. This 
method has recently been applied successfully to find the 
energy barriers for another organic molecular system^ 



TABLE I: Parameters for the atoms in the PTCDA system, 
as given in the AMBER force field 



atom 


AMBER atom type 


r* [nm] 


e[kcal/mol] 


C 


C/CA 


0.185 


0.12 








0.160 


0.20 





OS 


0.160 


0.15 


H 


H 


0.100 


0.02 



The calculation of the MEP for a PTCDA molecule 
from its minimum position within a PTCDA unit cell to 
the corresponding position in the neighboring unit cell 
reveals a diffusion barrier of Ed = 80 meV. The path 
and the potential energy along this path are shown in 
Fig. [2 Note that the minimum energy position of the 
molecule is not the same as the position in the bulk unit 
cell; this can be explained as the energy landscape for 
a single molecule should differ from that of a molecular 
layer. In a second calculation we determined the step 
edge barrier for a molecule. The model system was ex- 
tended by a fourth layer area on top of the three layers. 
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FIG. 2: Minimum energy path (MEP) for the center of mass 
(top) and potential energy along this path (bottom) for a 
PTCDA molecule on top of three layers of PTCDA in the 
(102) plane of the bulk a configuration. The molecules of the 
topmost layer are shown in grayshade. The MEP reveals a 
diffusion barrier for a single molecule of Ed = 80 meV. 



Then, again, the MEP for a PTCDA molecule over this 
single molecular step was calculated, revealing a fairly 
high barrier. Figure [3] shows the potential energy along 
the MEP over the step edge. The energy drops shortly 
before the step: The molecule gains energy when it is 
built in at the step. To go over the step, the molecule 
then has to overcome an energetic barrier of E$ = 750 
meV. 

To put these energy barriers into perspective, it is use- 
ful to introduce the concept of an onset temperature^ 
for a given kinetic process, which can be defined (some- 
what arbitrarily) as the temperature at which the pro- 
cess takes place at a rate of 1 s _1 . Assuming a standard 
attempt frequency^ 2 - of 10 12 s~ 1 , our estimated energy 
barriers imply an onset temperature of 33 K for diffusion 
and 315 K for step edge crossing. This indicates that 
the PTCDA molecules should be highly mobile at any 
temperature of interest, whereas interlayer transport is 
completely frozen: Even at 500 K, the ratio of the in- 
terlayer mobility to the intralayer diffusion rate is only 
e -(E s -E D )/k B T w 2x 10~ 7 , which implies that a molecule 
has to interrogate the step edge 5 x 10 6 times before de- 
scending. 

These predictions are consistent with growth ex- 
periments: When evaporated on an insulating KBr 
surface^ PTCDA exhibits three-dimensional growth. 
Even at low coverages, crystalline islands with a mini- 
mum height of 3 layers grow, and no wetting layers of 
PTCDA can be seen. On metal surfaces, the behavior 
is different!^ In this case Stranski-Krastanov-type be- 
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^ — ■ — minimum energy path (step crossing) 

| — ^— minim urn energy path (in- layer diffusion , see Fig. 2) 
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FIG. 3: (Color online) Potential energy along the MEP for a 
PTCDA molecule moving over a PTCDA step edge. The step 
edge barrier is Es = 750 meV. The MEP for in-layer diffu- 
sion is also drawn. Three important positions for a PTCDA 
molecule: (1) built into the step, (2) on the edge of the molec- 
ular step, (3) on the upper terrace. 

havior is found, in that the growth of three dimensional 
islands sets in after the formation of two complete lay- 
ers of molecules. The epitaxial structure of the first two 
layers, however, indicates a strong interaction with the 
substrate, which is not included in the calculations. The 
initial layer-by-layer growth does therefore not contradict 
our theoretical study. 

The small diffusion barrier we found also suggests the 
presence of mobile PTCDA molecules if experiments are 
performed at room temperature. Experimental evidence 
for this was found in scanning probe microscopy (SPM) 
experiments: When high resolution images of the top- 
most layer are taken, the images often become streaky. 



This is, e.g., described in ref4; the authors explain this 
behavior with mobile molecules that are diffusing on the 
topmost layer. 

Due to the high step edge barrier, molecules should 
be trapped on the topmost layer at room temperature. 
On the other hand, the small diffusion barrier might pre- 
vent molecules on top from nucleating and forming an 
incomplete additional layer. Indeed, in most SPM exper- 
iments, the top layers of PTCDA crystals seem complete 
at room temperature, as the mobile molecules cannot be 
resolved. In a recent experiment, 23 PTCDA molecules 
were found to form small clusters of molecules on top 
of crystallites when the sample is cooled down to 80 K, 
while at room temperature, the topmost layers seemed 
flat and complete. 

To summarize, we presented a study based on em- 
piricial potentials to determine the diffusion barrier Ed 
and the step edge barrier Es for an important organic 
molecule, PTCDA. On top of a cell of molecules in the 
bulk configuration, the MEPs for a single molecule be- 
tween minimum energy positions within neighboring unit 
cells and over a molecular step edge were calculated, 
yielding the estimates Ed = 80 meV and Es = 750 
meV. Although a very coarse model was set up, these 
results are consistent with experimental data. The accu- 
racy of our estimates of Ed and Es is difficult to assess. 
However, we note that the binding energy of a PTCDA 
molecule found using this approach, Eg 1 w 2.5 eV, is 
in close agreement with thermal desorption spectroscopy 
experiments^, which yield a binding energy of E e ^ p = 
2.2 eV. The method described in this paper should there- 
fore apply to many molecular systems and help to predict 
and explain effects in organic molecular crystal growth. 
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